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Although no experimental demonstration exists so far, sub-cycle electromagnetic pulses in the 
THz region, due to their time asymmetry, are considered to be one of the best options to orient 
molecules. However, we argue that orientation with half-cycle pulses, in the sudden regime, is 
an unphysical artifact that arises from a constant component of the field that is not a solution 
of Maxwell's equations. On the other hand, large instantaneous orientations can be achieved by 
time-symmetric single-cycle pulses of zero area. We show that, in this case, interferences among 
eigenstates of the time propagator give rise to a nonzero time average orientation, over the pulse 
duration, in spite of the existence of time-inversion parity. We also show that, due to the avoidance 
of the long-tail requirement, multi-single-cycle-pulse schemes can be designed to obtain a net average 
orientation over a rotational period. 

PACS numbers: 32.80.Lg, 33.80.-b, 42.50.Hz 



Rapid experimental developments of cooling methods 
and intense coherent sources has boosted an extraordi- 
nary interest in the understanding of phenomena related 
to the quantum control of molecular degrees of freedom 
[l|. In the last decade molecular alignment was one of 
the areas where impressive success was achieved, with 
the development of theoretical strategies and the experi- 
mental demostration of several techniques 0- However, 
molecular orientation has not progressed at the same 
speed. Control of molecular orientation with half-cycle 
pulses (HCPs) has been one of the most intensively stud- 
ied approaches in the last few years promising that high 
orientation can be obtained in a very efficient way [3|. 
Nonetheless, an experimental demonstration of the phe- 
nomenon has not been achieved yet. 

HCPs interact with the permanent dipole of molecules 
imparting them a sudden kick that transfers angular mo- 
mentum, and breaks rotational parity. A molecule in a 
rotational eigenstatc, is guided to a wave packet that 
contains rotational eigenstates of different parity, and 
for which the angle formed by the molecular axis and 
the field polarization vector is well defined. This view 
is based on a picture in which the HCP consists of a 
short and strong unipolar part during which a force is 
exerted on the molecule, and a weak long tail [J] which 
does not significantly alter the relative populations of the 
wave packet components. The long tail is required for 
the pulse to be propagable without considerable distor- 
tion because if the field area is not zero the pulse is not a 
solution of Maxwell's equations [f|. For a sudden inter- 
action the orientation achieved at the end of the unipolar 
part of the HCP is usually very modest, but during the 
long tail the phases of the components of the wave packet 
sinchronize producing at some later time a larger orien- 
tation. Schemes to maintain high molecular orientation 
have been devised based on sequences of pulses that im- 
part further kicks at each time the molecules start to lose 
their maximum orientation [6(. However, they only can 
work for heavy molecules for which the time to achieve 



the maximum orientation is comparable to the duration 
of the long HCP tail. Therefore, a new pulse can be sent 
at the correct time without interfering with the previ- 
ous one. On the other hand, light molecules go through 
several orientation/missorientation cycles during the tail, 
and before a new pulse kicks them the orientation is de- 
stroyed as shown in the upper panel of Fig. [T] A molecule 
with rotational constant B = 10 cm -1 loses its maximum 
orientation in less than 100 fs, while a molecule with 
B = 0.1 cm -1 remains oriented around 50 ps, which may 
be comparable to the tail duration. Thus, the long-tail 
requirement is an obstacle to obtain persistent orienta- 
tion except for heavy enough molecules. 

Apart from an apparently unnoticed paper by Sugny 
et al. Q , and a recent work by Fleischer et al. [|[ single- 
cycle pulses, that do not contain a long tail, have not 
been studied. Numerical evidence was presented in Q 
showing that symmetric pulses of zero area can produce 
significant instantaneous orientation during and after the 
pulse is over. In |8| it was experimentally demonstrated 
that intense single-cycle THz pulses induce field-free ori- 
entation that survives thermal averaging. In the present 
Communication I show that not only large peak orienta- 
tions but also a nonzero time average orientation during a 
single pulse of duration comparable to the rotational pe- 
riod can be obtained. This phenomenon is caused by the 
interference between pairs of eigenstates of the time prop- 
agator, as previously described for transport phenomena 
in quantum Hamiltonian ratchets 0, |T(| ■ Ratchets ex- 
hibit directed motion in the absence of an external bias. 
Recently, it has been shown that ratchet effects can oc- 
cur in coherent systems such as a Bose Einstein conden- 
sate realized with an optical lattice periodically modu- 
lated in time [TO]- Although this effect requires breaking 
all spatio-temporal symmetries simultaneously, transient 
currents can be generated in the presence of some sym- 
metry [l2j]. A similar mechanism is responsible for 
transient molecular orientation. 

Achievement of molecular orientation has been consid- 
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FIG. 1. (Color online). Instantaneous molecular orientation 
(cos 6) produced by THz pulses as a function of time, in re- 
duced units, for a molecule initially in the rotational eigen- 
state | J = 0,M = 0). Upper panel corresponds to half-cycle 
pulses with parameters a = 0.006(/i/B), T = 0M(h/B), 
ui = 200(B/h). In these units a rotational period is n. Solid 
line is for [iEq/B = 1000, and dashed line corresponds to 
fiEo/B = 100. When the unipolar part of the HCP ends no 
orientation is achieved. It is during the long tail that the 
molecules orient but taking into account that the average ori- 
entation, under field-free conditions, over a rotational period 
is zero, the orientation is lost if the tail lasts longer than the 
rotational period. In the lower panel the solid line shows ori- 
entation for the same initial state, for a zero-area pulse with 
HE /B = 12, a = 0A5{h/B), T = 3{h/B), u = 3{B/h). 
Dash-dotted lines give the orientation of the Floquet states 
that give a greater contribution to the initial state. 



ered traditionally more difficult |13j than the well un- 
derstood molecular alignment because it requires break- 
down of inversion symmetry. Three different schemes 
have been proposed, (i) combination of a strong nonreso- 
nant field and a weak static field [l4[ , (h) two-color laser 
schemes with resonant [ll| or nonresonant fill pu lses, 
and (iii) sub-cycle and few cycle THz pulses \vA Il8j . 
It is frequently argued that any orientational effect in 
molecules requires an external field asymmetric in time. 
Thus, the highly asymmetric temporal structure of HCPs 
is considered especially appropriate to achieve orienta- 
tion. Previous theoretical studies have been based on 
HCPs consisting of an unipolar part followed by an in- 
finite tail of infinitcsimally small negative intensity such 
that E(t)dt = 0. Since this tail does not interact 
with the molecule any possible orientation effect comes 
from the symmetric unipolar part, that can be written as 



E(t) — Eq exp(— (t — t ) 2 /a 2 ) cos(wi), where a is related 
to the Gaussian time width, and u) is in the THz region. 

Experimentally available HCP pulses of enough 
strength are short compared to typical rotational peri- 
ods. Thus, previous works were centered on the sud- 
den regime in which orientation is achieved only after 
the pulse is over. On the other hand, in the sudden 
regime, single-cycle pulses do not produce any notice- 
able orientation, but Sugny et al. pj showed that zero 
area pulses can produce large instantaneous orientations 
when its duration is not too short compared to the rota- 
tional period. Note, in passing, that longer pulses, in the 
adiabatic regime, can produce instantaneous orientation 
too. We show here that in the intermediate regime, but 
not in the adiabatic regime, the average time orientation 
during the pulse can be nonzero in spite of the presence 
of space-time symmetries. 

Constraints on the average orientation during a pulse 
or during several pulses that span a rotational period 
can be rigurousy understood by the existence of dynam- 
ical symmetries of the unitary time propagator. These 
symmetries can be analyzed more conveniently by con- 
sidering that the pulse is the first member of a periodic 
sequence of pulses, which allows us to study instead the 
symmetries of the Floquet Hamiltonian, as both opera- 
tors have the same eigenvectors [l9]. Thus, we will an- 
alyze the symmetries of U(T/2,—T/2) where T is the 
repetition period. The Floquet Hamiltonian (in reduced 
units [20]) for a rotating linear molecule in the presence 
of a periodic train of linearly polarized pulses, and period 
T is 
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where <p, the carrier envelope phase (CEP), is zero for 
single-cycle pulses of zero area, J is the angular momen- 
tum vector, B the rotational constant, /i the permanent 
dipole moment, and 9 is the angle between the polariza- 
tion vector of the field and the internuclear axis. Opera- 
tor J- acts in the extended Hilbert space where time, t' , 
is treated like another spatial variable [21| . Its eigenfunc- 
tions can be expanded in the basis of rotational functions 
(8, $| J, M) where M gives the projection of J along the 
polarization direction of the electric field, and $ is an- 
other Eulerian angle that defines the direction of J 



l^dj n exp{2mnt' /T) \ (9,§\J,M) . 

J>M V n / 

(2) 

Formally, eigenfunctions, x, of the Floquet operator U 
can be obtained from Eq. ([2]) by taking the projection 
t' = t = T/2. In practice, they can be easily calcu- 
lated using the generalized Floquet or (t,t') theory [22j], 
that gives the time evolution operator as the product 
]Im=o U i mT + T A (m - l)T/2] with 
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(3) 



where r is a small time interval, and F{&) is the matrix 
representing the Floquet Hamiltonian. For small enough 
t the matrix U can be truncated into a small matrix. 

Floquet theory gives us additional advantages over 
simpler propagation methods in terms of understanding. 
The method gives direct access to the matrix elements 
of the different field components, and the influence of 
a particular component of the electric field can be im- 
mediately understood. The Fourier expansion of a HCP 
pulse {(j) — 7r/2) contains a static component, which is 
the cause of molecular orientation for pulses in the sud- 
den regime as shown in Fig. [2j Dashed lines correspond 
to calculations in which the static component has been 
retained showing that large orientation can be obtained 
regardless of the pulse duration. However, real HCPs 
cannot have such a nonpropagable dc component. Far 
from the source, the dc and low frequency components 
are filtered out distorting the field [23|]. The solid lines 
in Fig. [5] are obtained after eliminating the effect of the 
static component, which is done by zeroing the terms 
with n — in Eq. ((3]). Panels a, and b in the fig- 
ure show that for short pulses the orientation is almost 
supressed. Thus, although previous theoretical studies 
had attributed molecular orientation by HCPs, in the 
sudden regime, to their asymmetric structure, the actual 
mechanism is the existence of an unphysical Stark effect. 
Longer HCP pulses have a dc component too, but it is rel- 
atively weaker than in the sudden regime. In this regime 
orientation is still possible even if the dc component is 
filtered out, as shown in panel c of Fig. [5J This com- 
ponent desymmetrizes the Floquet eigenstates, i.e., each 
one of the Floquet eigenstates that form the initial state 
may show a net average orientation during the pulse un- 
like the Floquet eigenstates for a zero-area pulse shown 
in the lower panel of Fig. [1] 

The electromagnetic field for a single-cycle pulse is 
an odd function of time, E(t) = —E(—t). Thus, the 
eigenvalue equation ,F(i')x(£') = is invariant un- 

der Stp ■ {x,9,t') -> (x*,7r - 9,-t'), which is the ro- 
tational version of the anti-linear transformation defined 
in The transformation 9 — > -k — 9, is a rotational- 

restricted version of the operation E* , that in molecules 
inverts the spatial coordinates of all the nuclei and elec- 
trons through the molecular center of mass j24(. Sym- 
metry Stp is called time-inversion parity [To| because it 
generalizes the notion of parity to the extended Hilbert 
space. For our system, it implies that if a given Flo- 
quet eigenstate is oriented (localized at 9q) at time —t 
it becomes antioriented (localized at 7r — 9q) at time t, 
resulting in a zero time average orientation during the 
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FIG. 2. (Color online). Orientation, (cos#), for HCP pulses 
(4> = 7r/2 in Eq. |T])) for three different regimes showing the 
effect of the static component of the field. Solid lines corre- 
spond to calculations in which the dc component of the field 
has been removed, while dashed lines have been calculated 
including it. The parameters, in reduced units, for thee three 
pulses are: (a) a = 0.006, T = 0.04, w = 200, fiE /B = 1000, 
(b) a = 0.07, T = 0.5, ui = 30, fiE /B = 100, and (c) 
a = 0.45, T = 3, cj = 3, fiE /B = 12. 



pulse. This dynamical symmetry does not constraint ei- 
ther the instantaneous orientations or the time average 
orientation during the pulse, except for Floquet eigen- 
states, contrarily to what has been implicitly assumed in 
previous studies. From a classical point of view, for each 
classical trajectory with orientation 9(t) there exists an- 
other one with orientation Tr — 9(t). Although the average 
orientation in the classical ergodic limit is zero Q , in the 
quantum case a nonzero average during the time span of 
one or several pulses can arise because nondiagonal ma- 
trix elements of the operator cos 6?, in the Floquet basis, 
can be nonzero. Thus, conservation of time-reversal par- 
ity does not forbid a net time average orientation if the 
initial state is a combination of several non-degenerate 
Floquet eigenstates as we show below. 

Taking as initial state a rotational eigenstate and ex- 
panding it in terms of Floquet eigenstates the time- 
evolved wave function can be written 



(t), (4) 



where b n = (xn(t = —T/2)\J,M), which gives for the 
instantaneous orientation 
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(^(t)\cose\^(t))=J2 ex M^n-e m )(t+T/2)]b* n b m O nm (t), 

nm 

(5) 

with O nm (t) = (Xn(t)\ COS 6\Xm(t)). 

Floquet eigenstates, in the extended Hilbert space, can 
be symmetric or antisymmetric with respect to time- 
reversal parity, x(^M) = = ' = X*( 7r — Using this 
fact and the property E*\J,M) = (-1) J+M \ J, M), it 
can be shown that the coefficients dj n in Eq. © are 
imaginary for J odd and real for J even or viceversa de- 
pending both on the symmetric or antisymmetric char- 
acter of the Floquet eigenstate and on the M value of 
the initial state. This implies for the coefficients of the 
eigenstates of U that c.j(t) = (-l) (J+fe) c}(-£), where 
c,j = '^2 n dj n ex.o(2itint/T). The integer k is zero for 
some eigenstates and one for others, resulting in O nn (t) = 
-O nn (-t), and O nm (t) = -(-l)^+ k ^O* m (-t). Tak- 
ing into account O nm {t) = O^ff), the time average over 
the pulse duration of (cos 9) (t) can be obtained integrat- 
ing exclusively over the second half of the pulse 

1 r T / 2 

((cos0)) = — / dt 2J A nm Q[O nm {t)exp[i(e n -e m )t], 

J° m>n 

(6) 

where 9[.] = Re[.] if k n ^ k m , and 9[.] = i x Im[.] if 
k n = fc,„, and the constants 

A nm = 2{6*6 m exp[i(e„ - e m )T/2] 



- (-l)^ +k ^b n b* m cM^m - e n )T/2}}. (7) 

There is no symmetry relation between (cos 9) (t) at 
any two different times greater than zero, and there- 
fore the integral in Eq. ^ can be nonzero in spite 
of time-reversal parity due to the interference between 
pairs of Floquet eigenstates encoded in the O nm terms. 
This mechanism can be loosely called transient breaking 
of time-reversal parity although, as we have shown, no 
symmetry is really violated. Fig. |3j shows the average 
orientation, during a rotational period, as a function of 
the field strength, for a sequence of ten pulses in panel a, 
and for a single pulse in panel b. For periodic sequences, 
equilibrium is reached after a few pulses since the average 
tends quickly to zero due to the growing phases in Eq. 

With current technology a strong single-cycle pulse can 
last no longer than a couple of ps, which seems to imply 
that obtention of a net average orientation over a time of 
the order of a rotational period would be possible by a 
single pulse only for the lightest molecules. However, the 
average orientation over a sequence of pulses is controlled 
by the differences between quasienergies unlike the aver- 
age orientation over a single pulse, which is controlled by 
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FIG. 3. (Color online). Time-average orientation, {{cos 6)), 
over a rotational period, for two different regimes, as a func- 
tion of field strength. Calculations were done by integration 
of Eq. ^ over the whole pulse. Parameters, in reduced units 
of time (h/B), are: (a) a = 0.045, T = O.Itt, and (b) a = 0.45, 
T = 7r. Thus, the time average in panel a corresponds to ten 
pulses and in panel b to one pulse. Panels c and d give the 
variation of the quasienergies with the field. Bigger (red) dots 
indicate the Floquet states with a greater contribution to the 
instantaneous wave packet. 



the Floquet eigenvectors. Thus, a nonzero time average 
can be achieved by using sequences that span a rotational 
period as shown in panel a of Fig. [3] The only condi- 
tion is being far from the sudden regime, otherwise the 
average orientation during each pulse will be very small. 
A greater stability in the average orientation is obtained 
with a sequence of pulses since the Floquet eigenvalues 
contributing to the instantaneous wave packets remain 
fairly constant with the field strength as shown in panel 
c of Fig. [3] On the other hand, for the longer pulse the 
existence of frequent avoided crossings (Fig. [3ji) changes 
often the character of the Floquet eigenstates that form 
the wavepacket, giving rise to larger variations on average 
orientation. Ben Haj-Yedder et al. i25j] explained that 
the most significant criterion for an efficiency-duration 
compromise regarding orientation phenomena, would be 
the maximization of the time average of (cos 9) over a 
rotational period. In the absence of an external field, 
orientation-dependent properties, averaged over the ro- 
tational period, vanish. Such efficiency criteria cannot be 
used for orientation strategies based on HCPs because in 
that case orientation is achieved under effective field-free 
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conditions, i.e., during the weak long tail. However, the 
criterion is meaningful for single-cycle pulses. 

Some confusion seems to exist in the bibliography con- 
cerning the symmetry properties that an external field 
must hold to orient a molecule. Two main consequences 
have resulted for the emerging field of THz quantum con- 
trol: (i) explanation of the orientation mechanism with 
HCPs based on their intrinsic time asymmetry, (ii) disre- 
gard for single-cycle pulses of ps duration and zero CEP 
as an option for orienting molecules. This is due, as we 
have shown here, to the fact that HCPs much shorter 
than the rotational period give rise to molecular orien- 
tation because they contain a strong static component 
that produce a Stark effect. This component is not a so- 
lution of Maxwell's equations and cannot propagate with 
the pulse, so the mechanism is an unphysical artifact. 
In the far field limit the pulse will be distorted and the 
dc component is filtered out. The resulting propagated 



pulse will then be unable to orient molecules. Contrar- 
ily, the regime where the pulse duration is comparable 
to the rotational period shows, for pulses of zero CEP, 
that the molecule gets sequentially oriented and antior- 
iented as the dipole force breaks parity. Interference 
among Floquet states causes a transient breakdown of 
time-inversion parity leading to significant average ori- 
entations for times similar or longer than a rotational 
period; a phenomenon impossible to obtain in the sud- 
den regime with half-cycle pulses. The average orienta- 
tion for a long sequence of pulses remains zero as corre- 
sponds to the equilibrium or long-time limit, due to the 
fact the interference terms accumulate large phases when 
time progresses 0. 
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